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ABSTRACT 


With the requirement for a spacecraft to maintain an orbital altitude band, a 
simple energy balance algorithm has been developed using a combination of radial 
distance and spacecraft specific energy for fixed-vector thruster control. While each 
trajectory produces a unique band, initial attempts at producing a pre-specified band have 
been unsuccessful. It is theorized that a certain radial bandwidth would correspond to a 
specific set of control parameters, and that by creating maps of the relationship between 
the two for various spacecraft configurations a method of maintaining the pre-specified 
band could be found. This thesis studies variations in spacecraft configurations and finds 
dependence of orbital bandwidth on thrust-to-drag ratio and ballistic coefficient. Also, 
within certain ranges of the control parameters, multiple trajectories produce equivalent 
radial bands. Analysis shows that all single-burn trajectories are characterized by similar 
efficiencies, and are less efficient than a Forced Keplerian Trajectory (FKT). 
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I. INTRODUCTION 


Optimizing orbit maintenance trajectories for low-Earth orbiting (LEO) spacecraft is 
increasingly important as space agencies around the world strive to maintain orbiting 
programs with shrinking budgets. Conceptually, an orbit may be maintained using a 
Forced Keplerian Trajectory (FKT), where thrust is used to cancel drag. This would prevent 
orbital decay, maintaining a true Keplerian orbit While this strategy is technically 
impractical due to limitations in thrust vectoring and thrust magnitude adjustments, it does 
provide a benchmark by which other strategies may be measured. While optimal control 
theory [Ref. 1], shows that an FKT is not the optimal solution to the orbit maintenance 
problem with respect to fuel, studies have yet to find a trajectory which is more fuel- 
efficient 

Historically, most orbit reboost strategies have been based on the Hohmann transfer, 
consisting of two thruster bums, where the first bum starts the spacecraft on an elliptical 
path to a higher altitude, and the second bum circularizes the spacecraft's orbit at the desired 
altitude. The Eccentricity Intercept Targeting and Guidance (EITAG) trajectory developed 
by Gottlieb [Ref. 2] notes variations of eccentricity, e, vs semi-major axis, a, and integrates 
backwards in time from the desired orbit A two bum trajectory is developed which reboosts 
and circularizes the orbit while minimizing errors in achieving the target orbit. 

A third orbit maintenance strategy proceeds from the notion that a spacecraft must be 
maintained within a specific radial band. A space station, for example, must be maintained 
within an altitude band that can be reached by various launch vehicles for servicing and 
manning requirements. A precise, circular orbit is not necessary. With this in mind, a single- 
bum trajectory, boosting the spacecraft as it nears the lower altitude boundary, then allowing 
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it to decay over time, would maintain the required band while reducing the number of 
required bums, and perhaps, the fuel required for orbit maintenance. 

Pauls [Ref. 3] and Wilsey [Ref. 4] developed a computer model for such a strategy, 
using various parameter variations to optimize the trajectory. While their results show that a 
single-burn trajectory is less fuel-efficient than the FKT benchmark, further study is 
warranted. In addition, an accurate method of predicting the actual controlled radial band is 
desired. 

By non-dimensionalizing the equations of motion, greater clarity and precision will be 
achieved when varying the control parameters. These results will be used to create maps of 
radius and energy which will allow better control of the spacecraft's radial band. 


2 




n. DEVELOPMENT OF THE ORBIT MAINTENANCE PROCEDURE 


A. NON-OPTIMALITY OF A FORCED KEPLERIAN TRAJECTORY 
While optimal control theory techniques show that a Forced Kepierian Trajectory 
(FKT) is not optimal with respect to fuel [Ref. 1 ], significant insight into the problem can be 
gained by looking into an extremely simplified version of the orbit maintenance problem 
[Ref. 5]. If an FKT is optimal, the theory should lead to an optimal control history where 
thrust is equal to drag. 



The problem may be simplified by stating that the only requirement of the control system is 
to maintain the orbit semi-major axis, a, and eccentricity, e, within such limits that the orbit 
remains within the specified band. These two orbital elements define the shape of the orbit 
within the orbital plane. This shape could theoretically vary from a circular orbit (tr=0) 
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where a varies from R m i n to Rm a « (see Figure 1), to aD elliptical orbit with perigee equal to 
Ratio and apogee equal to R MS , with e variation dependent upon the size of the band. 

Further simplification is possible by looking only at the portion of the problem relating 
to semi-major axis, a. This leaves a one degree-of-freedom system for optimization. To get 
an equation of motion for semi-major axis, the work-energy theorem is used 

AE-JFds (1) 


This can be used to relate die specific mechanical energy of the orbit to the forces involved, 
namely thrust and drag 

dE _ ( T-D ) ds 

dt ( m Jar 



Recalling the definition of the specific mechanical energy, e = -p / 2 a, 

de p da 
dt 2 a 2 A 


(3) 


This leads to an equation fora, 

*«2f!ll(T-D) 

pm 


( 4 ) 


The other equation of motion for the system comes from the rocket equation 


rh 


T 

Mo 


T 


V 

e 


(5) 
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where the specific impulse, Igp and the acceleration due to gravity are multiplied to define the 
exhaust velocity of the engine, v e . 

The final major simplification comes from focusing on the band to be maintained. 

When attempting to maintain a small enough band, the change in a is small, so the squared 
term is essentially constant. Additionally, such a band would require a small change in 
velocity, Av, to maintain it. Thus, v is also assumed to be constant. These assumptions allow 

Equation 4 to be rewritten as 

JL(T-D) (6) 

m 


where k is a constant. 

The optimal control problem is now formulated as finding the control history that best 
transfers the spacecraft from some initial condition a Q to some final condition a/. The 
optimality criterion is minimizing fuel burned, which can be stated as maximizing the final 
mass. The problem constraint is the available thrust, which is defined as 

0sT * T «, (7) 

Thrust is the single element of control in the simplified problem. From optimal control 
theory, the Hamiltonian is 


m v 

e 



(*) 
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where and \ m are the Lagrange multipliers, or costates of the system, corresponding to 
the semi-major axis and mass respeu ely. 

Pontryagin’s Maximum Principle (PMP) states that at every point on the trajectory of 
the system, the optimum control history maximizes the Hamiltonian. Stated mathematically, 

u’-argmax ufl ,H (9) 


By examining Equation 8, it is clear that the control variable, T, that maximizes the H- 
function depends upon the sign to the quantity in brackets, which is called the switching 
function, s 


m v 

e 


( 10 ) 


Noting the constraint upon thrust, application of the PMP shows that 

T. 


p- 

s>0 

(0 

s<0 


( 11 ) 


Called a bang-bang controller, this optimal control history contains no reference to drag; 
therefore the optimal control for such an orbital transfer does not define thrust equal to 
drag. The Forced Keplerian Trajectory is thus shown not be the optimal control method. 


B. ORBIT CONTROL STRATEGY 

It is desired to maintain a spacecraft within a specified radial band, using a single 
thruster firing to reboost the spacecraft when it nears the minimum radius. Wilsey [Ref. 4 ], 
shows that such a scheme is posable using radius and spacecraft specific energy as control 
parameters. Figure 1 shows such a radial band, defining radii and energies. 
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To facilitate this orbital control scheme, it is necessary to initiate thruster firing prior to 
the spacecraft decaying to the minimum altitude, and shut down thrusters prior to gaining 
the maximum altitude. Thruster firing is controlled by orbit radius, with 

6r-R„-r„ (12) 

where r t h is the radius of thruster initiation. Spacecraft specific energy is then monitored 
throughout the bum, with thruster cut-off occurring when 

E=E lh (13) 

with Etb being specific energy at thruster cut 'ft. By specifying values for 6r and En,, 
minimum and maximum radii will be achieved, denoted by r mD and r mx , respectively. This is 
pictured in Figure 2. 



The choices for ftr and F^, will produce a unique maximum band, but what values 
correspond to the desired band? Does a specific, desired band produce a specific trajectory? 
How do the bands and associated trajectories vary with changes in spacecraft parameters 
such asthrust-to-drag ratio, specific impulse, and ballistic coefficient? Ross, Pauls, and 
Wilsey (Ref. 6] propose a simple "energy-balance" algorithm where 6r and E| h are varied 
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until the desired radial band is achieved [Appendix. A |. By following this procedure and 
mapping 6r and E«h as they vary with the actual band attained for different spacecraft 
configurations, a relationship between the parameters is sought which would remove the 
need for iteration and allow orbit trajectory prediction. 

C. DEVELOPMENT OF THE EQUATIONS OF MOTION 

The initial investigation of a single burn, low-earth orbit trajectory will be simplified by 
assuming an orbit which is initially circular. The spacecraft is a non-lifting (blunt) body, 
therefore the only aerodynamic force which affects it is drag. The external forces thus 
perturbing the orbit are aerodynamic drag, gravity, and thrust. Looking at the orbital 
dynamics as a two-body problem, the coordinate system may be defined as shown in 
Figure 3. 



The equations of motion for this system are 



(14) 
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(15) 



where a r and a tr are the radial and transverse components of the inertial acceleration, XF r 
and 2F lr are the sums of the respective forces, and m is the spacecraft mass. As seen in 
Figure 3, drag and thrust can each be broken into components 

D r --Dsin(y) D r «-Dcos(Y) (16a&b) 

T r *Tsin(a) T tf -Tcos(a) (17a&b) 


where y denotes the flight path angle (between the velocity vector and the transverse axis), 
and a represents the thruster elevation angle (between the thrust vector and the transverse 
axis). Rewriting the equations of motion in polar coordinates using equations 16 and 17 
leads to 

2 u D T 

-L+—1 (18) 

r 2 m m 


0r+20f- 



m m 


(19) 


where p is the earth's gravitational parameter, and 8 is the angular position of the spacecraft. 

D. NONDIMENSIONALIZATION OF THE EQUATIONS OF MOTION 

/ 

When dealing with an orbital system, the use of standard dimensional units leads to 
large variations in numbers that cause computational errors. To counter this and to better 
study the effects of variation of certain parameters, a system of canonical units is developed 
to non-dimensionalize the equations of motion. 
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1. Definitions 

To nondimensionalize the equations, base units for length, mass, and time are 
defined. The base length is the radius of the target orbit, which is defined as the midpoint of 
the desired radial band, and is denoted by n>. The base mass is the initial spacecraft mass, 
and is denoted by mb- The base time is defined as the period of a circular orbit at a radius of 
it, divided by 2ji„ and is denoted by tb, 


t 


b 



( 20 ) 


These base units are then used to define the nondimensionalized variables for radius, mass, 
and time 


r 



( 21 ) 


m 



( 22 ) 


( 23 ) 

The variable 6 is dimensionless by definition. 

In addition to the three "basic” units, a base unit for density is also defined, and denoted 
by pb. Using a local exponential density model with atmospheric scale height $ 

( 24 ) 
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Pb is defined as equal to p G . Nondimensionalizing p with rb, and setting r re f equal to rb 
produces a nondimensionalized density 


p-e-w-o (25a) 

where 

P«K (25b) 


This is used in the determination of drag, as drag is given by 

D-i-pAC d v 2 (26) 

where A is the surface area, Cd the drag coefficient, and v the spacecraft velocity. A table of 
quantities, their nondimensionalized cousins and relating factors is given in Appendix B. 

2. Equation Nondimensionalization 

The above definitions are used to develop the nondimensionalized variable 
derivatives 

d 2 r _ d I dr \ dfdry^ _d_, p /b fl1 
dt 2 dt(dtj dt( d f/t7 dft b \ t 2 

dt 2 dt l dt J dt l df K dft/ \ t 2 


(27) 


(28) 


where primes denote a nondimensionalized parameter differentiated with respect to time. 
These derivatives are used to yield the nondimensionalized equations of motion 
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(29) 



r* 



rh m 


0T+20T 



m rh 


(30) 


where the nondimensionaJized values of both drag and thrust components are related to their 
dimensionalized counterparts by the same factor 

F - F , F - D,T (31) 

*Vb 


3. Parameter Development 

To study the effects of varying spacecraft design parameters on the orbit trajectory 
and their related bands, nondimensionalized parameters are developed which represent 
ballistic coefficient, engine specific impulse, and thrust-to-drag ratio. These parameters can 
be thought of as quantifying atmospheric effects, spacecraft shape, and engine sizing. 
a. Nondimensionalized Ballistic Coefficient, B 

Recalling equation 26, the nondimensionaliztion of drag requires the 
nondimensionalized density defined above, ^T, and a nondimensionalized velocity, v. 


Using the velocity equation 


v 2 »r 2 + r 2 0‘ 


and the relationships in equations 27 and 28 gives 


2 -2 ^”b 

V*-V _i 


(32) 


(33) 
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which is substituted, along with p into equation 26. Nondimensionaiizing drag using 
equation 31 gives 

D-1.J AC d v 2 £iIi (34) 

m b 

By defining the ballistic coefficient as the base mass divided by the product of the surface 
area and the drag coefficient, 


B 



(35a) 


the nondimensionalized ballistic coefficient is then defined as 

B —§_ 

P, r > 

which leads to the final nondimensionalized equation for drag 

D-j-JLil 

•* — 

• B 


(35b) 


(36) 


b. ManIThnut Parameter, p 

Fuel mass flow and thrust are related by the rocket equation, 

|f| __ T 

Mo 


(37) 


Following the example of equation 27 when nondimensionaiizing the derivative of mass and 
using the thrust relationship of equation 31, the nondimensionalized mass flow is written 
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ifi’.-T—— (38) 

By defining 

p- ~—- (39) 


variations in spacecraft specific impulse, 1^,, can be incorporated into a final 
nondimensionalized equation for both mass flow and thrust 

m' - - pf (40) 

c Thrust-to-Drag Ratio, 7\ 


The thrust-to-drag ratio is defined as the thrust divided by the drag at the target 
orbit. As both thrust and drag are forces and thus nondimensionalized by the same factors, 
this ratio may be written 



(41) 
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ID. COMPUTER MODEL DEVELOPMENT 


A. ORIGINAL DEVELOPMENT 

Originally developed by Pauls [Ref. 3] and refined by Wilsey [Ref. 4], a computer 
program, written in FORTRAN, was developed to simulate the orbital motion of a 
spacecraft. The main routine controlled input, output, and subroutine calls. Subroutines 
computed drag and orbital parameters, controlled thrust, and updated the equations of 
motion. A fourth-order Runge-Kutta numerical integration routine was used to integrate the 
equations of motion. The nondimen si on alized equations were simplified by defining new 
state variables, which transformed the two second order differential equations into a set of 
four first order differentials, 
a STATE VARIABLE DEFINITIONS 
A five element state vector, j, is defined as 

x,-r (42) 

(«3) 

x 3 -e (44) 

a de 

x 4 -e-~ (45) 

x 5 -tn (46) 
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These states can be nondimensionalized using the factors introduced previously and 
substituted into the equations of motion. Combining equations 29 and 30 with the 
nondimensionalized state variables produces 

dx, _ i d • T - 

— 7 T “ x i X 4 ~ t “ t" s,n Y +— s, na (47) 

at 5* m m 

i 



df 


2x 4 x, D f 

*■ - -—— cosy + -—— cosa 


x, x,m 


x, m 


(48) 


As the flight path angle, y, is a function of the velocity vector, the geometric functions can be 
represented by 

sin y-~ cosy-llfl (49a&b) 

v v 


Substituting these relationships into equations 47 and 48 leads to the final equations of 
motion as they are programmed into the simulation: 

d *. - 

-~-* 2 (50) 



% % 2 




mv 


+—sina 
m 


(51) 



(52) 



2x 4 x, Dx 4 f 


•cosa 


mv 


(53) 
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d *5 


dr 


-pt 


(54) 


C. PROGRAM EXECUTION 

Program execution begins with the input of variables from the user as we)) as an input 
file (see Appendix A). Following the definition of the initial conditions, the DRAG 
subroutine is called. The first call establishes the drag of the FKT benchmark. Execution 
then enters the main loop, where drag is called again, computing the instantaneous velocity 
as well as the density. With the velocity updated, the specific energy is computed and the 
CONTROL subroutine is called. This routine controls thruster activation and cut-off. With 
a thrust value thus determined, subroutine EQN updates the equations of motion. RK4 is 
the fourth order Runge-Kutta routine that integrates the state variables. Flight-path angle, 
fuel-burned totals and achieved altitudes are updated once every radian, while the routine 
ORBELMTS updates the classical orbital elements. The equations of motion are updated 
5000times per orbit, while all pertinent variables are sampled ten times per orbit. 


D. VALIDATION 

Since the program logic remains intact, the validation process is simplified. The initial 
baseline test consists of running the program with the spacecraft in an initially circular orbit 
with gravity the only external force. Radius, velocity, specific energy and angular 
momentum all remain constant, validating the basic equations without drag or thrust. 

The next step in the validation process compares the trajectory produced by the 
modified version with those produced by the earlier program for a consistent set of inputs. 
Since an in-depth validation is done by Wilsey fRef. 4], duplicating the results of the earlier 
version validates the modified program. The main difference between the two versions lay in 
the nondimensionalizing process, and the calculating of instantaneous spacecraft mass. To 
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achieve universality, all input variables and parameters are already nondimensionalized. 
Using input values which correspond to those used in the earlier study, the modified version 
produces equivalent results. 


E. DEVELOPMENT OF "SMART" VERSION 

The aim of this study is to determine values of 6r and E,h which will maintain a 

specified band, therefore a version of the program that accepts input values of 6r and Eih to 

simulate the trajectory is developed. This program is used as a test to validate the results 
obtained using the parameter mapping procedure. A listing of the program is found in 
Appendix C, while example input and output files are listed in Appendix O. 
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IV. PARAMETER MAPPING 

To remove the iteration on dr and Eu, which is required in the algorithm, a way of 

determining these control parameters for a given band prior to program execution is 
required. In addition, the variations in the control parameters for different spacecraft 
configurations must be known. 

To achieve these goals and possibly find a mathematical relationship between the 
parameters, a mapping technique is used. Running the simulation for a specific spacecraft 
configuration produces a set of trajectories. Multiple simulation runs for various 
configurations produces a family of trajectory curves. Combining all of the trajectories into 
a catalog where the achieved band is graphed against 6r and Eu, would bypass the 

algorithm's iteration process; for a given spacecraft configuration and a desired orbital band, 
values for 6r and Eu, could be pulled from the maps and input into the control program. 


A. PARAMETER VARIATION 

For thh work, a spacecraft configuration is defined as a set of the three 
nondimensionalized parameters previously introduced. Five values for each parameter are 
combined to produce a database of 125 different configurations. The parameter values and 
their dimensionalized counterparts are listed in Table 1. These values are chosen to 
represent a range of actual values for thrust-to-drag ratio, ballistic coefficient and specific 
impulse (ljp). Increasing values of T/D represent moving up to more powerful motors. 
Increases in specific impulse, the amount of thrust per unit weight of propellant, are 
reflected in decreasing values of p. The I$p values reflect those of different propellent 
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systems, from cold gas propellants (Ijp = 100), to ionic propulsion systems (Isp ~ 3000). 
Increasing values for B, the nondimensionalized ballistic coefficient, represent an increase in 


Table 1 PARAMETER VALUES FOR MAPPING PROCEDURE 



T/D 

B (B) 

P dsp) 1 

1 

50 

8.013e3 (50) 

7.902e0 (100) 

2 

75 

1,202e4 (75) 

3.161e0 (250) 

3 

100 

2.404e4 (150) 

1.581e0 (500) 

4 

250 

4.808e4 (300) 

7.902e-1 (1000) 

5 

500 

8.0l3e4 (500) 

2.634e-l (3000) 


spacecraft mass, decrease in surface area, or both, according to equations 35 a and b. A 

representative large orbiting platform with the following characteristics 

mass = 20,000 kg 
surface area = 60 m2 
drag coefficient = 2.2 

would have a ballistic coefficient of 150 kg/m?. 

a PROGRAM MODIFICATION 

Two modified versions of the computer simulation ate used to accomplish the mapping. 
One tracks the nondimensionalized orbital decay, 6r, while the other focuses on the thruster 

cut-off energy, E lh . 

1. Automatic variation of ordered band 

Both control parameters, 6r and E lb , can be defined by an "ordered" band. The 

amount of orbital decay defines the lower bound, while the thruster cut-off energy defines 
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the upper bound. By varying this ordered band, a range of values for fir and E lh can be 

simulated and mapped. A program loop added outside of all subroutine calls and data 
output commands automated repetition of the simulation, allowing the ordered band to vary 
from 50 km to 0 km (which would simulate the FKT), in 2.5 km increments. For each 
ordered band, the program simulated 50 orbits, or approximately three days of spacecraft 
flight 

2. "Bubble sort" of minimnm and maximum radii 

The actual band of a trajectory is defined as the radial distance between the 
minimum and maximum radii that is achieved during the simulation. Two 
nondimensionalized variables, and r^ are initially defined as equal to the initial orbit 

radius. As the program simulates the trajectory and the instantaneous orbit radius is 
computed, it is compared to the two variables. Whenever the instantaneous radius is less 
than the radius held in r Bm , the variable takes on this new, lesser radius. Conversely, r^ is 

updated to the instantaneous radius whenever the orbit goes above the current maximum. At 
the end of a trajectory simulation, f m and r K will hold the minimum and maximum radii 

achieved. The variables are reset for each trajectory, producing an actual band corresponding 
to each variation in fir and E (b 

3. Actual band vs. ordered band 

The purpose of the "ordered" band is tutorial; it provides a way to visualize and vary 
the control parameters fir and E |h . Once the simulation is run, the ordered band has no real 

significance. The actual radial band that is achieved for the given set of control parameters is 
ail that is of practical interest. An orbit maintenance problem as developed here is defined by 
the band to be maintained, not the individual obit radii. Any specified minimum and 
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maximum orbit radii form a radial band. By looking at the achieved band rather than the 
minimum and maximum radii separately, the dimension of the problem is reduced. 

C. MAPPING ACTUAL BAND VS. 6f AND EL 

IB 

Output from each of the parameter-mapping versions consists of the values of the 
respective control parameter and the corresponding actual band produced. An example plot 
for thruster cut-off energy is shown in Figure 4. This depiction of the trajectory is then 
combined with similar trajectory curves corresponding to different spacecraft 
configurations, producing a catalog of trajectories. With a specific spacecraft and a required 
band to be maintained, values of the control parameters 6r and E th can be pulled from the 

plots and used to maintain the spacecraft's orbit, removing the need for iteration in the 
proposed algorithm. 



Figure 4 Example map of achieved band vs. thruster 
cut-off energy 
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V. ANALYSIS AND RESULTS 


Analysis of the parameter maps consists of investigating the effects of spacecraft 

parameter variation on achieved band. While 125 different spacecraft configurations are 

catalogued, investigation is focused on four configurations which best exhibit the noted 

phenomena. The following configurations are studied: 

Config 533 Config 234 Config 342 Config 453 

T/D = 500 T/D = 75 T/D=100 T/D = 250 

B = 150 kg/m 2 B=150kg/m2 B=300kg/m2 B = 500kg/m2 

Isp = 500 sec Isp = 1000 sec 1^ = 250 sec Isp = 500 sec 

The prefix D- or E- on the maps indicate the source file of the data; D represents actual 

band plotted against orbital decay {6r), while E represents actual band plotted against 

thruster cut-off energy (E lh ). The first digit refers to T/D ratio, the second to B, and the third 


to Isp. On each plot, only one parameter is varied. 

Accuracy of the mapping technique is also tested by using the control parameters in the 
"smart” version of the simulation and noting the size of the actual band achieved. All figures 
discussed are located in Appendix E. 


A. VARIATION OF SPECIFIC IMPULSE, Lp 

Figures 1 through 4 show the effects of varying a spacecraft engine's specific impulse. 
In all cases, the effect of varying Isp is extremely small. Figure 1 shows no variation between 
trajectories. It is also noted that the curves for actual band vs. 6r exactly match the shape of 
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the actual band vs. E u . This is due to the way the control parameters are varied together, 

using the concept of the "ordered" band introduced in the previous chapter. 

Figures 2,3, and 4 all show families of trajectory curves corresponding to spacecraft 
whose configurations differ only in the value of their engine's specific impulse; thrust-to- 
drag ratio and ballistic coefficient remain constant within each figure. While individual 
trajectories differ, the general shape in each case remains the same. Also, as the actual band 
narrows, the trajectories become practically identical. This convergence represents identical 
control parameters producing equivalent bands for different spacecraft configurations. 

Recalling the definition of the mass/thrust parameter p (equation 39) and Table 1 it is 
clear that despite large changes in specific impulse, the parameter p remains a small number 
that varies little. Figures 1 through 4 confirm that variations in specific impulse have little 
effect on the orbital band maintained using the single-burn algorithm. 

B. VARIATION OF BALLISTIC COEFFICIENT 

Figures 5 through 8 show the effects of varying the spacecraft's ballistic coefficient. 
Again, the values erf thnist-to-drag ratio and (in these cases) specific impulse remain 
constant throughout each figure. Comparing the trajectories in any of the individual maps 
shows considerably more effect on actual band by variation erf ballistic coefficient than was 
produced by varying specific impulse. 

Generally, the size erf the actual band for a given control parameter decreases with 
increasing values erf the nondimensionalized ballistic coefficient. Figure 5 shows this trend 
with only small deviations. Figures 6 and 7 show larger variations between individual 
trajectories, bid the shrinking of the actual band does continue. Figure 8 shows less 
variation, which points towards a relationship between the ballistic coefficient and the thrust- 
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to-drag ratio. Varying ballistic coefficient has less of an influence on configurations with 
larger values of T/D. 

Recalling equations 35 and 36, increasing spacecraft mass or reducing spacecraft 
surface area leads to an increase in ballistic coefficient, which in turn decreases drag. As the 
thrusters modeled in this simulation operate at a constant thrust, those configurations with 
large thruster values feel little effect from the decrease in drag, while those with smaller 
thrusters, and thus smaller T/D ratios, more easily "sense" the decreased drag. With thrust 
and drag smaller, the algorithm is better able to control the actual band. 

C. VARIATION OF THRUST-TO-DRAG RATIO 

Figures 9 through 12 show the effects of varying the spacecraft’s thrust-to-drag ratio. 
Each plot represents a family of spacecraft with consistent values of ballistic coefficient and 
specific impulse, but different T/D ratios. Comparing individual trajectory curves within a 
map shows that increasing thruster size generally increases the actual orbital band achieved. 

As with ballistic coefficient variation, the effects of varying T/D are more profound than 
when changing specific impulse. Once again, the effect of large thrusters is most noticeable, 
as they seem to cancel out irregularities present in configurations possessing smaller 
thrusters. In each figure, the trajectory corresponding to that configuration with the largest 
thrust-to-drag ratio (T/D - 500) is a smooth curve, and in some cases almost a straight line. 
Larger values for T/D ratio inhibit the single-burn algorithm; for a given value of either 
control parameter, a larger thruster produces a larger band. Smaller thrusters, which may be 
thought of as being more responsive to this control scheme, produce a tighter orbital band. 

D. MULTIPLE TRAJECTORIES FOR A SPECIFIED BAND 

An assumption made during the formulation of the single-bum algorithm and the 
mapping procedure covered here deals with the relationship between the control parameters 
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and the band they produce. Given a spacecraft and a set, periodic thruster firing schedule 
(defined values of 6r and E lh ), a certain radial band will be maintained. The parameter 

mapping technique essentially reverses this process, starting with the required radial band 
and looking at the required thrusting scheme. 

It seems reasonable to assume that since a unique thruster firing schedule produces a 
specific radial band, the converse will also hold; a certain radial band will yield unique 
values of 6r and E th . As can be seen in Figures 2 through 12, this is not the case. On each 

map, there is a region where a specific band can be produced by multiple values of each 
control variable. Such a choice between control parameter values could lead to new 
optimizing opportunities: which trajectory is "best*? To investigate this development, the 
"smart" version of the simulation is used to study different thruster firing schedules for a 
specific spacecraft configuration. 

E. "SMART" PROGRAM VALIDATION 

Figure 13 shows die two {dots relating actual band to the control parameters for 
Configuration 234introduced previously. For nondimensionalized radial bands between 
approximately .0076and .0095 times the base orbit radius, there are three possible values 
for both control parameters. To validate the smart program and investigate the control 
parameter relationships, the following trajectories are chosen which produce nearly equal 
radial bands: 

Trajectory Actual band 
Traj.l 0.006610134 

Traj.2 0.006606818 

Traj3 0.006789484 





The values of 6r and E lk corresponding to each trajectory are used by the smart program to 


generate the specific simulation for 50 orbits. 

1. Trajectory comparison 

The orbit radius of each individual trajectory is shown in Figures 14 through 16. A 
combined plot is shown in Figure 17. While each trajectory produces a radial band of 
comparable width, the bands and the thruster firing patterns are quite different Figure 14 
shows Trajectory 1 's decay to initial thruster firing taking just over ten orbits, followed by a 
long bum. After thruster cut-off the orbit's eccentricity is quite significant, as seen in the 
rhythmic variation in nondimensionalized radius (Rbar) as the orbit decays again to r th . 


Figure 15 shows a similar pattern for Trajectory 2, with a smaller orbital decay and a shorter 
bum-time. 

Figure 16 shows a much different pattern. Decay to the initial thruster firing occurs 
rather quickly, with the orbit proceeding directly to a highly eccentric state. The eccentricity 
of the orbit causes the spacecraft to hit the specified thrusting radius as well as the thruster 
cut-off energy more frequently. The orbit actually becomes an ellipse whose apogee and 
perigee are bounded by the minimum and maximum radii of the orbital band. Noting the 
minimum and maximum radius values for Trajectory 3 plotted in Figure 16 and the 
definition of eccentricity 




r - r 


r +r 


(55) 


the eccentricity is found to be approximately 0.0043. Figure 18 compares the eccentricities 
of the three trajectories, and that of Trajectory 3 confirms the computed result, as it tends 
towards a steady-state value of approximately 0.0045. 
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2. Efficiency comparison 

Given these differing bum patterns, which one is best? A performance index must 
be defined to be optimized; this study looks at maximizing fuel efficiency. Figure 19 shows 
the mass of fuel burned for each trajectory, as well as that for a Forced Keplerian 
Trajectory. 

It is clear that the conclusions made by Wilsey [Ref. 4] are confirmed; none of the 
single-bum trajectories are more fuel-efficient than an FKT. What is of interest, however, is 
that the average slope of the fuel-burned curve appears independent of thruster bum pattern. 
Regardless of the frequency and length of thruster firing, all three single-bum algorithm 
produce the same average fuel bum rate. 

3. Variation Of Orbital Simulation Length 

In comparing these trajectories, it is helpful to look at different orbital parameters, 
including spacecraft specific energy. The nondimensionalized specific energy is defined as 



and is shown for each orbit in Figure 20. The near vertical increases in specific energy 
correspond to thruster firings; the greater the rise in energy, the longer the bum. Trajectory 
3 quickly reaches a state where the specific energy is tightly controlled by a high-frequency, 
short-bum schedule. With only two bums showing for the other trajectories, however, 
sound conclusions cannot be drawn; more data is required. An additional series of 
simulations is run for a period of 100 orbits, with orbit radius, specific energy, and fuel 
mass burned shown in Figures 21 through 24, respectively. 


28 



a. Transition to a Tight Energy State 

From Figure 21, it is seen that Trajectory 2 undergoes a significant change around 
its 65th orbit; the radial hand maintained expands significantly. Comparing the three 
trajectories in Figure 22, it is seen that Trajectory 2 transitions to a tightly governed energy 
state very similar to that reached earlier by Trajectory 3 (although the energy level of 
Trajectory 2 is higher). This transition appears due to the eccentricity of the mbit increasing 
to a point where the orbit grazes the prescribed band; the perigee and the apogee of the orbit 
are consistently in contact with die bounds defined by the control parameters, as shown in 
Figure 23. Despite the change in the thrusting pattern and the increase in radial band size. 
Figure 24 shows that the slope of the fuel-burned curve remains consistent, above that of the 
FKT. 






VL CONCLUSIONS AND RECOMMENDATIONS 


The purpose of this thesis was to create a catalog of control parameter maps for various 
spacecraft configurations, using a proposed single-bum trajectory algorithm for low-earth 
orbit maintenance. Such maps are used to remove the iteration required in the algorithm 
(Appendix A), assuring maintenance of the desired band and removing the need for human 
intervention in the control process. Additionally, study of the relationships between control 
and spacecraft configuration parameters led to new insights into the orbit maintenance 
problem. Finally, comparison of all single-bum trajectories to the Forced Keplerian 
Trajectory (FKT) benchmark for optimization of fuel-efficiency could not corroborate the 
optimal control theory finding that the FKT is not the optimal trajectory. 

The FORTRAN simulation developed by Pauls [Ref. 3] and refined by Wilsey [Ref. 4] 
was further modified by nondimensionalizing the entire simulation, then developing two 
versions to accomplish the parameter mapping. A "smart" version was also developed which 
would accept as input the control parameters mapped out in the thesis, thus testing the 
validity of the (dots. 

Three nondimensionalized spacecraft configuration parameters related to thrust-to-drag 
ratio, ballistic coefficient, and specific impulse were defined and varied through a range 
practical values. Their individual effects upon the spacecraft's trajectory were studied. 
Variations in specific impulse produced little change in the trajectories. Variations in 
ballistic coefficient and thrust-to-drag ratio produced more pronounced changes between 
trajectories, with a coupling noted between the two. While increasing the value of the 
ballistic coefficient generally decreased the size of the maintained band, the influence of the 
ballistic coefficient was significantly reduced for high thrust values. High values of T/D 
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produced larger radial bands, and also seemed to cancel out any orbital irregularities present 
in spacecraft with smaller, less powerful thrusters. 

An interesting discovery was that a specific radial band could be produced by different 
combinations of control parameters; the relationship between an orbital band and a thruster 
firing schedule was not one-to-one. This allowed for further optimization studies, as a 
certain single-bum scheme might actually be superior to an FKT. 

While the F1CT was shown to be consistently better than all single-bum trajectories in 
this study, the differences in thruster firing patterns proved intriguing. Different trajectories 
producing comparable radial bands had drastically different thruster firing times and 
frequencies. High orbit eccentricity, where the orbit's perigee and apogee were bounded by 
the radial band, sent the trajectory into a periodic, high-frequency, low bum-time thrusting 
scheme, where specific energy varied little. Fuel-bum curves for each single-bum trajectory 
maintained a consistent slope, even after the orbit transitioned to this tightly governed 
specific energy state. 

This thesis opens the door for further study as well, since some questions were 
unanswered, and new questions have arisen. A comparison of the single-bum algorithm to a 
dual-bum algorithm such as EITAG [Ref. 2] is warranted. While not superior to an FKT, a 
single-bum orbit maintenance algorithm may still be better than the dual-bum solution. 
Additionally, what effects do other orbital parameters such as flight-path angle and 
spacecraft velocity have upon the trajectory? Is the transition to the tight energy state 
globally valid? And on a grander scale, what is the optimal orbit maintenance strategy? What 
form of bang-bang control will minimize fuel burned and provide the most efficient means 
of maintaining a continued presence on orbit? With man's continued presence in space 
seemingly contingent upon maximizing performance and efficiency while minimizing cost, 
such study is certainly warranted. 
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APPENDIX B 
NondimensionaJization Table 


Quantity (=) 

Parameter (*) 

Factor 

radius : r 

RBAR 

rb 

mass : m 

MBAR 

mb 

time : t 

TBAR 

t b 

velocity : v 

VBAR 

rb/tb 

energy : E 

EBAR 

p/n> 

thrust : T 

THBAR 

mbTb/tb 2 

density : p 

RHOBAR 

Pb 

drag : D 

DBAR 

mbTb/tb 2 

ballistic coeff : B 

BBAR 

PbTb 

mass flow : MDOT 

MDOTBAR 

mb/tb 

atm scale ht : p 

BETABAR 

1/n, 

angular momentum : h 

HBAR 

n» 2 /tb 

semi-mqjor axis : a 

ABAR 

rb 

thrust/drag parameter : Td 

TDBAR 

TBAR/DOBAR 
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APPENDIX C 


A. PROGRAM LISTING 


PROGRAM ORB MAI NT 

ORBMAINT is a self-contained, nondimensionalized orbit maintenance program 
that will control the orbit of one body orbiting another. The program consists of a main 
routine controlling input, output, and subroutine calls, and five subroutines which 
compute atmospheric drag, the equations of motion, the standard orbital elements, and 
the number and duration of thruster firings. A 4th order Runge-Kutta integrator is used 
to solve the equations of motion. The program also performs comparative calculations 
for a Forced Kepierin Trajectory (FKT) 


VARIABLE DEFINITIONS 
ARRAYS: 

XBARO) 

XBAR(2) 

XBAR(3) 

XBAR(4) 

XBAR(5) 

XBDOTO) 

XBD0T(2) 

XBD0T(3) 

XBD0T(4) 

XBD0T(5) 

RLE INPUT: 

ROBAR 

MOBAR 

FORBIT 

STEP 

PRNT 

CETABAR 

TDBAR 

BBAR 

PBAR 

USER INPUT: 

ALPHA 

ALPHAR 

DLTRBR 

ETHBR 

COUNTERS: 

TBAR 

INDEX 

KOUNT 


nondimensionalized orbit radius 
nondimensionalized orbital radial velocity 
theta 

nondimensionalized orbital angular velocity 

nondimensionalized spacecraft mass 

derivative ofXBAR( 1) 

derivative of XBAR(2) 

derivative of XBAR(3) 

derivative of XBAR(4) 

derivative of XBAR(5) 

initial nondimensionalized orbit radius 
initial nondimensionalized spacecradft mas 
number of orbits to simulate 
time increment step size 
print interval step size 
atmospheric scale height 
nondimensionalized thnist-to-drag ratio 
nondimensionalized ballistic coefficient 
nondimensionalized mass/thrust parameter 


thrust angle 

nondimensionalized orbital decay: ROBAR-RTHBR 
nondimensionalized spacecraft energy 

current time (initial value=0) 

(initial value=0) 

(initial value=0) 
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INITIAL AND BOUNDARY CONDITIONS: 

DOBAR initial nondimensionalized atmospheric drag 

EOBAR initial nondi mensionalized spacecraft energy 

VOBAR initial nondimensionalizedspacecraft velocity 

RMNBR nondimensionalized minimum orbit radius 

RMXBR nondimensionalized maximum orbit radius 

PROGRAM VARIABLES: 


ABAR 

APOGEE 

DBAR 

EBAR 

ECC 

GAMMAD 

GAMMAR 

HBAR 

FUELKEP 

TOTKEP 

MBRFUEL 

PERIGEE 

PRDBAR 

RHOBAR 

RTHBR 

TBARF 

THBAR 

THBRCAP 

THBRKEP 

THGMA 

VBAR 

CONSTANTS: 

PI 

J 


nondimensionalized semi-major axis 
nondimensionalized orbit apogee 
nondimensionalized atmospheric drag force 
nondimensionalized spacecraft energy 
orbital eccentricity 

flight path angle 

nondimensionalized angular momentum 
keplerian fuel burned per program iteration 
total keplerian fuel burned 
nondimensionalized mass of fuel burned 
nondimensionalized orbit perigee 
nondimensionalized orbit period 
nondimensionalized atmospheric density 
nondimensionalized radius of thruster firing 
simulation stop time 
nondimensionalized spacecraft thrust 
nondimensionalized spacecraft thrust capability 
nondimensionalized keplerian thrust 
flight path angle at thruster initiation 
nondimensionalized spacecraft velocity 

P* 

state variable index 


variable declarations: 

IMPLICIT REAL*8(A-H,M-Z) 
DIMENSION XBAR(5),XBDOT(5) 


constant definitions: 

PI=DATAN( 1.0D+00)*4.0D+00 
J=5 


initialize keplerian fuel totals 
FUELKEP=0.0D+00 
TOTKEP=0.0D+00 
CCCCCCCCCCCCCCCCCCCCCC 
C 

OPEN( 10, FILE='input', STATUS='OLD') 
OPEN(20, FILE='orbpar\ STATUS='NEW’) 
OPEN(30, FILE='orbel\ STATUS= , NEW) 
OPEN(40, FILE=’orbetc\ STATUS='NEW’) 
C 
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C read input file 

REACH 10,1) R0BAR,M0BAR,FORBIT,STEP,PRNT,BETABAR, 

*TDBAR,BBAR,PBAR 
1 FORMAT* 10(/,21X,D 15.9)) 


convert time and intervals to radians 
TBARF=2.0*P1 *FORBIT 
STEP=2.0*PI*STEP 
PRNT=2.0* PI * PRNT 


read user inputs of thrust angle and control parameters 
PRINT*, 'ENTER ALPHA' 

READ*, ALPHA 
ALPHAR=ALPHA*PI/180.0D+0 

PRINT*, 'ENTER DLTRBR’ 

READ*, DLTRBR 

PRINT*, 'ENTER ETHBR' 

READ*, ETHBR 


initialize counters 
TBAR=0.0D+00 
INDEX =0 
KOUNT=l 


initialize state variables 
XBAR(1)=R0BAR 
XBAR(2)=0D+0 
XBAR(3)=0D+0 
XBAR(4M 1/R0BAR)**1.5 
XBAR(5)=M0BAR 


initially, set achieved orbit radii to center of orbital band 
RMNBR=R0BAR 
RMXBR=ROBAR 


set radius for thruster initiation 
RTHBR=R0BAR-DLTRBR 

CALCULATIONS 


call DRAG to determine initial and keplerian atmospheric drag 
CALLDRAG(BBA1LXBAR3ETABAR,VBAR,DBAR) 

DOBAR=DBAR 

C commence calculation loop, call DRAG to determine s/c velocity and atm drag 
100 CALLDRAG(BBAR,XBAR3ETABAR,VBARJ)BAR) 

C determine s/c energy 

EBAR=(VBAR*VBAR)/2.0-( 1/XBAR( 1)) 

C determine spacecraft thrust capability 
THBRCAP=TDBAR*D0BAR 
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C determine kepierian thrust (=atmospheric drag) 

THBRKEP=DOBAR 

C call CONTROL to fire or shutdown thrusters 

CALLCONTROL(ETHBR,XBAR,THBARJRTHBR,EBAR,THBRCAP) 

C call EQN to update the orbital equations of motion 

CALLEQN(DBAR^CBAR v XBDOT,VBAR,THBARALPHAR,PBAR) 

C call RK4 to do the runge-kutta integration over time 

CALL RK4(TBAR,XBAR,XBDOT J ,STEP,INDEX) 

check for completion of runge-kutta integrator, if not done, loop back to 100 
IF (INDEX .NE 0) GO TO 100 

runge-kutta complete: update variables for computation of orbital elements 
and output 

GAMMAR=ATAN(XBAR(2)/(XBAR( 1 )*XBAR(4))) 
GAMMAD=GAMMAR* 180.0D+0/PI 
compute flight path angle during thrusting, THGMA 
IF (THBAR .GT. 0.0D+00) THEN 
THGMA =GAMMAD 

ELSE 

THGMA=0.0D+00 
END1F 

RBAR=XBAR(1) 

MBAR=XBAR(5) 

MBRFUEL=M0BAR-MBAR 
FUELKEP=THBRKEP*PBAR*STEP 
TOTKEP=TOTKEP+FUELKEP 
IF (RBAR .LE RMNBR) RMNBR=RBAR 
IF (RBAR .GE RMXBR) RMXBR=RBAR 
ACTBAND=RMXBR-RMNBR 

call ORBELMTS to determine the classical orbital elements 

CALLORBELMTS(EBARBBAR,VBAR,GAMMAR,ECCABARAPOGEE 
♦PERIGEEPRDBAR) 

check if print interval acheived. if not, skip ouput section 
IF (KOUNT .LT. DNINT(PRNT/STEP)) GO TO 200 

OUTPUT 

WRTTE(2D,2)TBAR/(2.0*PI),RBAR,VBAR,MBRFUEUEBAR 

2 FORMAT(2XE6.2^XEl 1.9,2X,F11.93X.F11.9,2X,F11.9) 

C 

WRITE(303)TBARy(2.0*PI)ABARECCAPOGEEBERIGEEPRDBAR 

3 FORMAT(2XE62,lXEl 1.93X JB63XE11 -9,1X,F11,90X,F8.4) 

C 

WRITE(40,4)TBAR/(2.0*Pi)T)BAR,GAMMAD,THBAR,TOTKEP 

4 FORMAT(1XB633XE12.9,1XE10.4,1XE11.9,1X,15.12) 

C 

C once printing completed, reset counter 
KOUNT=0 
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C update counter 

200 KOUNT=KOUNT+1 

C 

C check if simulation complete, if not loop back to 100 
IF (TBAR .LT. TBARF) GO TO 100 
C 

PRINT*, 'Actual band achieved: ' 

PRINT*, ACTBAND 
C 

END 

CCCCCCCCCCCCCCCCCCCCCC 

C DRAG calculates the spacecraft velocity and the atmospheric drag acting on the 
C spacecraft, using a local exponential atmospheric density model 
C 

SUBROUTINE DRAG<BBAR^CBAR,BET ABAR, VBARJDBAR) 

IMPLICIT REAL*8(A-H,M-Z) 

DIMENSION XBAR(5) 

C 

C calculate spacecraft velocity 

VBAR=<(XBAR(2)*XBAR(2))+(XBAR( 1 )*XBAR(4))**2)**0.5 
C calculate nondimensionalized atmospheric density 
RHOBAR=EXP(-BETABAR*(XBAR( 1)1)) 

C calculate nondimensionalized atmospheric drag 

DBAR=RHOBAR*VBAR*VBAR/(2.0*BBAR) 

C 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCC 

C EQN updates the spacecraft equations of motion 

C 

SUBROUTINE EQN(DBAR,XBAR,XB DOT, VBAR,THBAILALPHAR,PBAR) 

IMPLICIT REAL*8(A-HM-Z) 

DIMENSION XBAR(5) ,XBDOT(5) 

C 

PI=DATAN( 1 .OD+00)*4.0D+00 
C 

A=XBAR( 1 )*XBAR(4)*XBAR(4) 

B= 1 /(XBAR( 1 )*XB AR( 1)) 

C=DBAR*XBAR(2V(XBAR(5)*VBAR) 

E=THBAR*S1N(ALPHAR)/XBAR(5) 

C 

F=2.0*XBAR(4)*XBAR(2VXBAR( 1) 
G=DBAR*XBAR(4)/(XBAR(5)*VBAR) 

H=THBAR*COS(ALPHA R)/( XB A R( 5) * XB AR( 1)) 

C 

XBDOT( 1)=XBAR(2) 

XBDOT(2)=A-B-C+E 

XBDOT(3)=XBAR(4) 

XBDOT(4)=-F-G+H 

XBDOT(5)=-PBAR*THBAR 

C 
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RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCC 

C ORBELMTS calculates the nondimensionalized orbital elements 
C 

SUBROUTINE ORBELMTS(EBAR,RBAR,VBAR,GAMMARJ5CC,ABAR, 
*APOGEEJ > ERlGEE,PRDBAR) 

IMPLICIT REAL*8(A-H,M-Z) 

C 

PI=DATAN( 1.OD+00)*4.0D+00 
C 

C calculate angular momentum 

HBAR=RBAR*VBAR*COS(GAMMAR) 

C calculate eccentricity 

PROB=( 1.0+2.0*EBAR*HBAR*HBAR) 

IF (DABS(PROB) .LT. 1.0D-12) THEN 
ECC=0.0D+00 

ELSE 

ECC=PROB**0.5 

ENDIF 

C calculate semi-major axis 

ABAR=-1.Q/(2.0*EBAR) 

C calculate apogee 

APOGEE=ABAR*( 1.0+ECC) 

C calculate perigee 

PERIGEE=ABAR*( 1 .(LECC) 

C calculate period 

PRDBAR=2.0*PI*( ABAR** 1.5) 

C 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCC 

C RK4 does the fourth order runge-kutta integration 

C 

SUBROUTINE RK4(TBAILXBAR*XBD0T r IJSTEP,INDEX) 

IMPLICIT REAL*8(A-H,M-Z) 

INTEGER INDEX.I 

DIMENSION XBAR(5)*XBDOT(5),SAVED(5),SAVEX(5) 

C 

INDEX=INDEX+1 
GO TO (1,23,4)4NDEX 
C 

1 DO 101=14 
SAVEX(I)=XBAR(I) 

SAVED(I)=XBDOT(l) 

10 XBAR(I)=SAVEX(IH0.5EM)*STEP*XBDOT(I) 

TBAR=TBAR+0.5D+0*STEP 
RETURN 
C 

2 DO 20 I=1J 

SAVED(I)=SAVED(I)+2.0D+0*XBDOT(I) 
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20 XBAR(I)=SAVEX(I)+0.5D+0*STEP*XBDOT(I) 

RETURN 

C 

3 DO 301=1 J 

SAVED(I)=SAVED(I)+2.0D+0*XBDOT(I) 

30 XBAR(I)=SAVEX(I)+STEP*XBDOT(I) 

TBAR=TBAR+0.5D+00*STEP 

RETURN 

C 

4 DO 401=1 J 

40 XBAR(I)=SAVEX(I)+STEP/6.0D+0*(SAVED(I)+XBDOT(I)) 

INDEX =0 
RETURN 
C 

END 

CCCCCCCCCCCCCCCCCCCCCC 

C CONTROL determines if the spacecraft is within radial and energy limits prescribed 
C by the desired radial band and controls firing of the thrusters accordingly 
C 

SUBROUTINE CONTROUEITffiR^GB AR,THBARJITHBRJEBAR, 
THBRCAP) 

IMPLICIT REAL*8(A-H,M-Z) 

DIMENSION XBAR(5),XBDOT (5) 

C 

C check if thrusters already on 

IF (THBAR .EQ. THBRCAP) GO TO 100 
C check if orbit radius is less than or equal to thrusting radius 
IF (XBAR( 1) .LE. RTHBR) THEN 
C check if s/c energy is less than or equal to thrusting energy 

IF (EBAR .LE ETHBR) THEN 
C if both are true, then fire thrusters 

THBAR=THBRCAP 

ELSE 

C if not, do not fire thrusters 

THBAR=0D+0 

ENDIF 

ELSE 

THBAR=0D+0 

ENDIF 

C 

C check if s/c energy is greater than thrusting energy, if so, turn off thrusters 
100 IF (EBAR .GE ETHBR) THBAR=0D+0 
C 

RETURN 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
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APPENDIX D 


A. SAMPLE INPUT FILE 


)D+00 
)D+00 
)D+02 
)D-04 
)D-01 
I.407286740D+02 
7.500000000D+01 
2.403901513D+04 
7.902309754D-01 


l.( 

1.C 

1. ( 

2 . < 
l.< 


iiiiiiii 




niiiitt 


hi hi 11 




ROBAR 

MOBAR 

FORBIT 

STEP 

PRNT 

BETABAR 

TDBAR 

BBAR 

PBAR 


B. SAMPLE OUTPUT FILES (through 20 orbits) 


SPACECRAFT CONFIGURATION: 234 
T/D:75 B: 150 Isp:1000 

REQUIRED BAND: 0.006035 

DLTRBR: 0.001035681 

ETHBR: -0.499482695 

ACHIEVED BAND: 0.0059946 


FILE: ORBPAR 


P RB .rr. 

1.00 

2.00 

3.00 

4.00 

5.00 

6.00 

7.00 

8.00 

9.00 

10.00 

11.00 

12.00 

13.00 

14.00 

15.00 

16.00 

17.00 

18.00 

19.00 

20.00 




0.999733807 

0.999457389 

0.999169910 

1.003238879 

1.003357471 

1.003534507 


1.003680096 


1.003794299 
1.003874691 
1.003919779 
1.003930433 
1.003904591 
1.003843996 
1.003747894 
1.003617000 
1.003452807 
1.003257664 


QtUt 


1.002785952 


1.002513219 


VELOCITY 

1.000133063 

1.000271285 

1.000415031 

0.997268162 

0.997127460 

0.996949890 

0.996803384 

0.996688313 

0.996606844 

0.996560448 

0.996548591 

0.9965728S9 

0.996631972 

0.996726219 

0.996855001 

0.997016952 

0.997209856 

0.997430684 

0.997677081 

0.997947130 


fuelbi 


Ulttiat 




0. 

0. 

0. 

0.002423820 


M.l 11 M 111 


I.UIIUM 


-300133172 

-300271584 

-300415663 


-.499499684 


0.002644306 

0.002913851 

0.003182619 

0.003451645 

0.003720414 

0.003988924 

0.004257950 

0.004526460 

0.004795487 

0.005063997 

0.005332249 


rjV5T» 


0.005869011 

0.006138037 

0.006406547 

0.006674541 


-.499522178 

-.499523400 

-.499524904 

-.499526248 

-.499527664 

-.499529162 

-.499530408 

-.499531863 

-.499533080 

-.499534523 

-.499536089 

-.499537673 

-.499539165 

-.499540455 

-.499543844 

-.499543844 



RLE- ORBEL 


..— \ 


1.00 0.999733728 0.000002 0.999736079 0.999731376 6.2807 

2.00 0.999457127 0.000005 0.999461975 0.999452278 6.2781 

3.00 0.999169364 0.000008 0.999176871 0.999161858 6.2754 

4.00 1.000956558 0.003012 1.004017094 0.997986175 6.2926 

5.00 1.000956558 0.003013 1.003972874 0.997940243 6.2922 

6.00 1.000954109 0.003008 1.003965144 0.997943075 6.2922 

7.00 1.000951095 0.003003 1.003957068 0.997945122 6.2922 

8.00 1.000948404 0.002999 1.003950002 0.997946805 6.2921 

9.00 1.000945566 0.002995 1.003943187 0.997947945 6.2921 

10.00 1.000942563 0.002991 1.003936143 0.997948982 6.2921 

11.00 1.000940067 0.002988 1.003930476 0.997949658 6.2920 

12.00 1.000937151 0.002985 1.003924458 0.997949844 6.2920 

13.00 1.000934713 0.002982 1.003919724 0.997949702 6.2920 

14.00 1.000931822 0.002980 1.003914742 0.997948902 6.2920 

15.00 1.000928684 0.002978 1.003909285 0.997948082 6.2919 

16.00 1.000925510 0.002975 1.003903572 0.997947449 6.2919 

17.00 1.000922520 0.002973 1.003898246 0.997946794 6.2919 

18.00 1.000919936 0.002972 1.003894303 0.997945570 6.2919 

19.00 1.000916822 0.002971 1.003890236 0.997943408 6.2918 

20.00 1.000913145 0.002969 1.003884904 0.997941385 6.2918 


RLE- ORBETC 
ORBIT DRAG 


1.00 0.000021599 
2.00 0.000022462 
3.00 0.000023396 
4.00 0.000013114 
5.00 0.000012893 
6.00 0.000012571 
7.00 0.000012313 
8.00 0.000012114 
9.00 0.000011975 
10.00 0.000011899 
11.00 0.000011880 
1X00 0.000011924 
13.00 0.000012028 
14.00 0.000012194 
15.00 0.000012424 
16.00 0.000012718 
17.00 0.000013078 
18.00 0.000013500 
19.00 0.000013968 
20.00 0.000014543 


GAM MA T HRUST 


- 0.0001 
-0.0003 
-0.0004 0. 

0.1157 0. 

0.1045 
0.0688 
0.0721 
0.0546 0. 

0.0365 0. 

0.0179 
-0.0009 0. 

-0.0197 0. 

-0.0382 
-0.0563 
-0.0737 
-0.0902 0. 

-0.1056 0. 

-0.1197 0. 

-0.1324 0. 

-0.1434 0. 


JK TE UEL 



0.000103273109 
0.000206546218 
0.000309819327 
0.000413092436 
0.000516365544 
0.000619638653 
0.000722911762 
0.000626184871 
0.000929457980 
0.001032731089 
0.001136004198 
0.001239277307 
0.001342550415 
0.001445823524 
0.001549096633 
0.001652369742 
0.001755642851 
0.001858915960 
0.001962189069 
0.002065462178 





Band vs. Delta Rbar 
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Band vs. Ethbr 
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Band vs. Delta Rbar 



Band 







Band vs. Ethbr 



Band 












Delta Rbar 

















Band vs. Ethbr 
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Figure 3b 
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Band vs. Ethbr 
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Band vs. Delta Rbar 
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Figure 13b 
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